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ABSTRACT 

We have generated complementary halo mass estimates for all groups in the Galaxy 
And Mass Assembly Galaxy Group Catalogue (GAMA G'^Cvl) using a modified caus- 
tic mass estimation algorithm, originally developed by |Diaferio &: Geller ( |1997 ). We 
calibrate the algorithm by applying it on a series of 9 GAMA mock galaxy light cones 
and investigate the effects of using different definitions for group centre and size. We 
select the set of parameters that provide median-unbiased mass estimates when tested 
on mocks, and generate mass estimates for the real group catalogue. We find that 
on average, the caustic mass estimates agree with dynamical mass estimates within 
a factor of 2 in 90.8 ± 6.1 % groups and compares equally well to velocity dispersion 
based mass estimates for both high and low multiplicity groups over the full range of 
masses probed by the G'^Cvl. 



1 INTRODUCTION 



A quantitative understanding of the largest structures in 
the Universe provides a rigorous test of cosmology and dark 
matter simulations (Eke et al. 1996V One of the principal 



members, e.g., Hughes (19891; Carlberg et al. (19961; Gi- 



rardi et al. (19981; Tucker et al. (20001. By assuming that 



assumptions of the ACDM paradigm is that galaxy groups 
serve as tracers of the underlying dark matter haloes. It is 
thought that prior to the decoupling of matter and radia- 
tion, dark matter particles, with their smaller interaction 
cross-section, formed overdense regions under the influence 
of gravity. Following decoupling, baryonic matter is left to 
free fall into these overdense regions and form the building 
blocks of the large scale structure of the universe ( |SpringeI| 
et al.|2005 1 . Empirical measurements of the distribution and 
meisses of galaxy groups therefore provide a powerful con- 
straint for different dark matter models. 

Traditionally, mass estimates of groups are calculated 
virially, via measurements of the velocity dispersions of their 



the group is in virial equilibrium, the dynamical mass of the 
group follows the relation M cx a^R. The obvious limita- 
tion of this method is that the relation will hold only out 
to the virial radius of the galaxy group, so mass estimates 
made using galaxies beyond this radius become less reliable. 
More accurate mass estimates at large radii are possible us- 
ing weak lensing (Kaiser et al. 19951, however the obvious 



drawback of this approach is the observational challenge in- 
volved in measuring the lensing signal for a large galaxy 
group sample. 

A different approach to group mass estimation is to look 
at the distribution of galaxies within a group in redshift 
space (the projected distance r from the group centre and 
the line-of-sight velocity v with respect to the median group 
redshift for every member of the group) and estimate the 
group escape velocity by interpreting the distinct shape of 
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this distribution. This method of analysing galaxy group 
members in redshift space was first introduced as a mass 



estimator by Diaferio fc Geller (19971 and Diaferio ( 1999 1 



(hereafter DG97 and D99) and is known as the caustic mass 
estimation technique. An early version of the method was 
first used in an attempt to constrain cosmological parame- 
ters, specifically Q.m ( Rcgos fc Gcller_^1989j with little suc- 
cess. In identifying that velocity measurements of galaxy 
groups were heavily affected by random motions, as well as 
by comparing observed caustics to those predicted by cosmo- 
logical models, DG97 determined that this method is unre- 
liable for determining Q.m- Navarro et al. (19971 went on to 



demonstrate that the density profiles of dark matter haloes 
do not vary with cosmological parameters. However, the am- 
plitude Air) = min{|u„|, ji;;!} of these caustics, defined as 
being half the difference between the upper and the lower 
line-of-sight velocities Vu and vi does provide a measure of 
the gravitational potential 4>{r) of the group. The method 
has been subsequently used by Diaferio et al. (20051; Diafe- 



(20091; Serra et al. (2011 1 as a robust way of calculating 



masses of galaxy groups using all of the positional informa- 
tion from a group (out to -R200) and as a test for cluster 
membership ( Rines fc Diaferlo| [2006). Additional research 
has shown that for a set of rich X-ray luminous clusters, 
caustic mass estimates agree to within a ratio of 1.03 ± 0.11 



with mass estimates obtained via lensing analysis (Rines 
|et al.|2003| |. 

In this work we apply the caustic method to the Galaxy 
and Mass Assem bly (GAMAQ Galaxy Group catalogue 
(hereafter G^Cvl; Robotham et al. 20111 and attempt to 
provide complementary caustic mass estimates to the dy- 
namical mass estimates of the group haloes within the cata- 
logue with the aim of verifying these masses. By also apply- 
ing the method to mock light cones that mimic the GAMA 
data we are able to carefully calibrate our algorithm to 
produce median-unbiased mass estimators for each galaxy 
group using only redshift and positional information out to 
radii that are well beyond the virial radius. 

The GAMA project ( [Driver et al.|2009[[20Tl] ) is an on- 
going major galaxy survey covering 21 bands of the elec- 
tromagnetic spectrum from the ultraviolet through to radio 
wavelengths using a multitude of ground and space based 
telescopes. The aims of the survey are to cover a region of 
~ 360 deg'^ and to obtain ~ 400,000 redshifts for galaxies to 
a magnitude of tab = 19.8 mag. During phase one of GAMA 
we observed three fields of 4 x 12 degrees centred at a = 9h, 
5 = 1 deg (G09), a = 12h, 5 = deg (G12) and a = 14.5h, 5 
= deg (G15). In phase two we are expanding these to 5 x 12 
degrees as well as gathering data in three new fields. One 
of the principal scientific goals of GAMA is to understand 
and better constrain the halo mass function (HMF 
fc Schechtef] [19741 |Lacey 



Cole||1994i [Sheth et al.] 
Warren et al.|2006[ ' |T'inker et al.|2008| ), which describes the 



Press 
200TT 



relationship between the mass and number density of dark 
matter haloes. To that end, it is very important for any 
catalogue of galaxy groups produced to have accurate and 
reliable mass estimates, particularly for low-mass groups. 

This paper is structured as follows: in Section 2 we 
briefly describe the theoretical framework behind the caus- 



http://www.gama-survcy.org/ 



tic mass estimation method, our changes to the algorithm, 
and a step-by-step description of its implementation. In Sec- 
tion 3 discuss the data set used in this work (mock galaxy 
group catalogues created for GAMA and the actual groups 
themselves), and our results are presented in Section 4. We 
summarise our results in Section 5. Throughout this paper, 
consistent with the cosmology used to create the GAMA 



mocks and Robotham et al. (20111, we use a cosmology of 
= 0.25, = 0.75, Ho = h lOOkms"^ Mpc'^ 



2 METHODOLOGY 

The caustic method relies on analysing the distribution of 
group members in a redshift space diagram, which is defined 
in D99 as the plane (r, v) of the galaxies, where r is the pro- 
jected radial separation from the group centre, and v is the 
line-of-sight velocity relative to the group centre of mass. 



The spherical infall model of Regos & Geller ( 1989 1 predicts 



the existence of two trumpet shaped 'lines' on this plane 
where the phase-space density in redshift space is infinite; 
and in practice these trumpets are observed when looking at 
both simulated and real groups. By definition, galaxies out- 
side of these caustics are beyond the turnaround radius of 
the group. If we assume that galaxies lying outside the caus- 
tics are considered to be escaping the group, it follows that 
the caustic describes the escape velocity vl^^ir) of a cluster 
as a function of distance r from its centre. The term 'caustic' 
refers to the formation of a singularity at a location where 
the Jacobian of a co-ordinate transformation vanishes; in the 
case of the caustic method, this is the transformation from 
real to redshift space. In the spherical model, this causes the 
galaxies within a group to collapse to the peak redshift of 
the group along the line-of-sight, placing them somewhere 
between the centre of the group and the turnaround radius. 
In redshift space this translates to the trumpet-like lines 
that describe the escape velocity of the group. 

Here follows a brief review of the physical justification 
behind the caustic method. Full details of the model can 
be found in DG97, D99 and ( [Serra et ar][2011| ). Assuming 
a spherically symmetric model, the escape velocity within 
a shell of radius r for a group is given by v'^{r) = — 20(r). 
Given our position as observers, it is the line-of-sight compo- 
nent (uios) of this escape velocity that determines the loca- 
tion of the caustics, but this value depends upon the escape 
velocity profile of the cluster which we may not always know. 
Instead, we require an expression for the caustic amplitude 
that is independent of the escape velocity profile. [Serra et ah] 
( 2011 1 determine such an expression via the thought process 



summarised below: 

We begin by looking at the velocity anisotropy parame- 
ter I3{r) — 1 — {{v'g) + {v^^))/2{vl) where ve, and Vr corre- 
spond to the longitudinal, azimuthal and radial components 
of an individual galaxy's velocity respectively. If we assume 
that cluster rotation is negligible, («|) = {v^) = (wfos) and 
(Wr) = {v^) — 2{vi^^) where v\as is the line-of-sight com- 
ponent of the velocity. Rearranging this for {v^) and in- 
corporating the equality of all velocity components into the 
expression for /3(r) gives: 



3 - 2/3(r) 
l-/3(r) 



= ("L)ff('-) 



(1) 
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where 



9{r) = 



3 - 2/?(r) ^ 2{«L) + (vi) 
1 - P{r) 



(2) 



Given that the potential of a system is related to its 



escape velocity in the form of —2(j> - 



-(r)), it is possible 



to link the potential to the caustic amplitude if we make the 
assumption that A'^{r) — (w^sc los) in the form 



- 2<^(r) = {v^sc,ios) 9{r) = A {r)g{r) 



(3) 



which reduces our problem to having only one unknown: the 
parameter /3. The final link in the chain is to consider the 
mass of an infinitesimal shell: 



Gdm : 



-2(t){r)T{r) dr = A' {r)g{l3)F{r) dr 



(4) 



with J-(r) — 



-2TTGp(r)r-' 
4>(r) 



Combining the expression for J-{r) 



and G dm brings us back to the familiar result for the mass 
of a shell; Gdm — 4-KGp{r)r^ dr. Integrating Eq. |4]gives 



(5) 



GM{<r)= / A\r)g{r)F{r)dr 



This expression is close to what we need, but is limited 
by the fact that the density profile of the system needs to 
be known in order to get a value for the mass. To overcome 
this, D99 assume that J-fiir) = J^(r)g{j3) is a slowly varying 
function with respect to the radius of the system in hierar- 
chical clustering scenarios, and this result is confirmed by 



Serra et al. (20111. We therefore set J-pir) to be constant 
and adopt the value of — 0.7 of |Serra et al.| ( |2011[ ) (this 
is not critical, as we later adjust our mass estimates by a 
scaling factor, discussed below), giving the final expression 

GM{<r)=Tis f A^{r)dr (6) 
Jo 

which can be used to provide mass estimates for our data 
sample. 



2.1 Caustic mass estimation algorithm 

A successful caustic mass estimation algorithm must be able 
to successfully infer the continuous dark matter distribution 
in the halo from a discrete set of points determined by galax- 
ies. The most important goal of our algorithm is to correctly 
determine the location of the caustic for a group. To do this, 
we project the group into an area of redshift space and gen- 
erate a kernel that describes the continuous density of the 
group within this area ( |Pisani|1993| ). Based on this density 
distribution the algorithm then determines a threshold at 
which the caustic is placed. 

For a given galaxy group for which member positions 
are known (a, S, and z), the projection into redshift space 
takes place via the following transformations: 

cDa{zc) , , 

r= — - tanip (7) 



and 



Ho 



1+Zc 



(8) 



where Da is the comoving distance to the galaxy, Zc is the 
redshift of the group centre, and -i/; is the angular separa- 
tion of a member galaxy from the group centre at redshift 



z along the line of sight. The area of redshift space is there- 
fore determined by the full radial extent of the group and 
its range of line-of-sight velocities. Consider A'' galaxies in 
a cluster distributed in a redshift diagram with coordinates 
X = (r, v). Using an adaptive kernel method (Silverman 



19861, we describe the density distribution of these galaxies 



/,(x) 



JV 

-T- 

N ^ 

1=1 



where K is the adaptive kernel 



K(t) = 







ty iit<i, 

otherwise 



(9) 



(10) 



and hi = hchoptXi is the local smoothing parameter. \i = 
[tZ/iC'^O]^^^ where /i is equation where he = \i = 1 
for any i and log 7 = log [/i(xi)] jN . The motivation for 
using an adaptive kernel estimator is to have a density es- 
timator that can adapt to density distributions where the 



true probability density changes quickly ( Pisani 1993 1; this 



is generally a caveat of fixed kernel estimators that risk to 
oversmooth or undersmooth the probability distribution. Fi- 
nally, the optimal smoothing parameter ftopt is 



3.12 

7V1/6 



„2 I „2 



1/2 



(11) 



where and (j„ are respectively the uncertainties in the 
galaxy coordinates. The positional uncertainty Or is calcu- 
lated from the astrometric uncertainties in GAMA and is 
negligible, while for we use the same value of 55 km sT^ 



as used in Robotham et al. (2011 1 



Performing this calculation can take a great deal of 
time, particularly when it comes to calculating jq and opti- 
mizing for the best value oih^. A faster, time saving way of 
obtaining the density estimator is to use Fast Fourier Trans- 
forms to convolve a two-dimensional histogram of the data 
in redshift space with the adaptive kernel ( Silverman|1986 1 . 
Other alternative implementations of the algorithm include 
using friends-of- friends algorithms and binary trees, but in 
our case, where the galaxy groups have already been selected 
and statistically well-defined, the FFT approach works best. 
The process is described in |Silverman| ( [1986) for one dimen- 
sion, but can easily be extended to two. 

We begin by creating a two-dimensional normalised his- 
togram of the galaxies in redshift space (projected radius 
from the group centre versus line-of-sight velocity), and over 
the same parameter space, a histogram of calculated values 
for the kernel given in Eq. ( |10[ |. The smoothing parameters 
Ai, ^opt and he are used to adjust the values of the data 
histogram. The density estimate /^(r, v) is defined as the 
inverse FFT of the product of the forward FFT of the data 
and kernel histograms, or in other words: 



/q(r, u) = JT-i (jr(data) x J'(kernel)) 



(12) 



where T and denote forward and inverse fast Fourier 
transforms respectively. The end result of this calculation is 
a two-dimensional matrix that describes the density distri- 
bution of the group galaxies in redshift space; note that we 
use the modulus of this matrix in order to discard any phase 
shifts caused by the Fourier transforms. The matrix has di- 
mensions of 2® X 2*; we find that any size lower than 2"^ x 2^ 
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does not provide our algorithm enough resolution to give 
reliable results. The caustics are drawn on the density dis- 
tribution on locations where fq{r, v) — k, and k is obtained 
by minimising the function S{k, R) taken from D99: 



S{K,R) = \{vL).,R-4{v^)Rf 



(13) 



where the term {vcsc)k,r is the square of the mean escape 
velocity at R for a given value of k. This corresponds to 
the average size of the caustic amplitude from the group 
centre to the maximum projected radial distance R for a 
given value of k and {v'^)r is the group velocity dispersion 
taken from the G^Cvl. To minimise the function S{ii,R), 
we wish to find the value of k for which the average caustic 
size is equal to 4{i;^)_b. We use the R function OPTIM to 
do this; OPTIM is a general-purpose optimisation function 



based around the Nelder-Mead algorithm (Nelder & Mead 



19651, which lends itself particularly well to this task as 



S{k, R) is a parabolic function with only one well-defined 
minimum. Once the location of the contour is determined 
(black line in Figure [l| , the algorithm draws the caustics 
(green lines in Figure [Tralong min{|?;„|, \ vi\} where and Vi 
correspond to the upper and lower values of the line-of-sight 
velocity of the group along the contour. The algorithm scans 
through the density distribution in bins of r, and for each 
bin selects the minimum of these two velocities {vu{r) and 
vi{r), and reflects it along the line-of- sight velocity axis. The 
caustic amplitude beyond the maximum extent of the group 
is artificially set to 0, even though often the caustic closes 
before the maximum radial extent of the group is reached 
(see Figure [T|. 

Based on mock catalogues of galaxies built with N-body 
simulations in DG97, there is a constraint on the logarithmic 
derivative of the caustic amplitude: dln^/dlni? ^ 2. Any 
values of A{r) for which this derivative does not hold are 
considered to be the result of the caustic algorithm coming 
up with the wrong location for the caustic at that partic- 
ular radius, often due to excessive foreground/background 
galaxies. Instead, in these cases we use a value for the caustic 
amplitude such that dln^/dln_R 
(20TT1). 



1/4, 



Serra et al. 



To summarise, our algorithm works as follows: 

(i) Convert the galaxy positions in redshift space into a 
two-dimensional histogram of galaxy number densities, 

(ii) create another histogram of the same dimensions con- 
taining values for the kernel as per Eq. ([9}, 

(iii) calculate fq{r, v) using Eq. (121, 

(iv) calculate the best value for k with Eq. (131, 



(v) fit the caustics by reading off the minimum value of 
\vu\ and \vi\ along r, whilst ensuring that the derivative in- 
equality holds, 

(vi) integrate between the caustics to estimate the mass 
of the group using Eq. ( |14[ | and scale it accordingly to obtain 
a median-unbiased estimator. 

It is important to stress that this FFT implementation of the 
kernel estimator still retains the calculations and expressions 
given in this section; we simply speed up the procedure of 
calculating the density distribution by applying the kernel 
with a FFT. 

Our final mass expression is 





T — I — I — I — I — I — I — I — I — I — r- 



—I — I — I — I — I — I — 



I I I I I I I I I I I I I I I I 1 I I I I I I I I I I I I I I I I I I 11 I I 



I I I I I I I I I I I I I I I I 
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Projected distance (h Mpc) 



Figure 1. Four examples of the placement of the contour k 
(black) and the caustics fitted to it (green) for four FoF mock 
galaxy groups (whose galaxies are shown as the black points) 
from the G^Cvl in descending order of luminosity. It is evident 
here that the FFT method used to estimate the density distribu- 
tion causes the final caustics to be very smooth with respect to 
the caustics drawn in DG97 and D99. This introduces a source of 
error in our caustic mass estimates. 



0.7A, 



/i-iMq G/(M0km2s- 



^Mpc) 



Airfdr 



(14) 



where r is given in units of Mpc and A in s~^ km. 

A^{r)dr is calculated by discretising A{r) over a set of 
equally spaced steps and Ac is the caustic mass scaling fac- 
tor. In Figure[l]we show example caustic fits for four friends- 
of-friends mock galaxy groups of descending total luminosity 
(from 10^^ h-2 Lq to 10^ h'^ Lq). 

Despite the computational efficiency of the FFT method 
to calculate the density estimate, there are a number of 
drawbacks that need to be adressed. To begin with, the area 
over which the 2D histograms for the data and the kernel 
are created need to be larger than the area the data spans. 
This is to avoid the kernel (and the resulting density es- 
timate) wrapping around the borders due to the periodic 
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nature of Fourier transforms and results in a very smooth 
density distribution compared to those shown in DG97 and 
D99. This will effectively artificially increase the sizes of our 
caustics and cause a systematic overestimation of the caus- 
tic mass (though the scaling factor Ac corrects for this); the 
smoothing is independent of group size as shown in Figure 
[l] We do not consider the presence of background galaxies 
in this analysis, and make the assumption that the friends- 



of-friends algorithm used in Robotham et al. (20111 has re- 



covered group members as accurately as possible. This as- 
sumption is tested in Section 3. 



3 DATA 

3.1 The GAMA Galaxy Group Catalogue 

The groups presented in the G^Cvl have been identified 
using a slightly modified friends-of-friends algorithm ( |Press| 
fc Davis]|1982 Huchra fc Geller||1982 I to correctly account 
for the distortions due to peculiar velocities. It achieves this 
by considering projected separations independently of radial 
positions. Figure [2] shows cone diagrams containing galaxies 
and groups for the G12 field, where it is possible to see how 
they trace the large-scale-structure. 

The final group catalogue contains approximately 
14,000 galaxy groups with 2 members out to tab = 19.4 
(fields G09 and G15) and 19.8 mag (G12 field) which encom- 
pass about 40% of the galaxies in the full GAMA catalogue. 
Robotham et al. estimated their masses by using the rela- 
tion M Qc a^R where a and R correspond to the velocity 
dispersion and the projected group radius. The group ve- 
locity dispersion is calculated using the GAPPER algorithm 



given in Beers et al. ( 1990 1 where all the recession velocities 



for A'^ galaxies within a group are ordered, and the velocity 
dispersion is estimated using calculated values for the gaps 
between velocity pairs Vi+i — Vi for i — 1, 2..., A^ — 1. Group 
radius estimates were calculated for radii that contain 50%, 
68% and 100% of the galaxies within each group (Radso, 
Radgg and Radioo respectively), and the projected group 
centre was located by three different methods, discussed be- 
low. The mass of each group was estimated by 



A 



2 

"FoF 



/i-iMq G/(M0km2s-2Mpc) ^km^s-^ / ft-^Mpc 



Radj 



(15) 



where G is the gravitational constant (6.673 x 10~^" km^ s~^ 
Mpc), RadpoF and (Jfof are the radius and velocity disper- 
sion of the group obtained by the methods described above. 
A is a scaling factor that is required to obtain a median- 
unbiased estimate of the friends-of-friends mass with respect 
to the real halo mass. It varies depending on multiplicity 
and median redshift of a group. In this work we mimic the 



multiplicity and redshift subsets of Robotham et al. (20111 



to ensure a direct comparison. We consider every galaxy in 
a group when running our algorithm, while the dynamical 
mass estimates do not contain the full radial galaxy density 
profile information. The projected group centre is defined 
in three different ways in the G^Cvl: the first centre corre- 
sponds to the rAB luminosity centre of light (CoL) of all the 
galaxies associated with the group. The iterative group cen- 
tre is estimated by calculating the tab CoL of the cluster, 
and then rejecting the most distant galaxy from this centre. 



This process is then repeated until only two galaxies remain, 
at which point the brightest one of these is selected as the 
group centre. In the final method, the brightest group galaxy 
in the rAB-band is selected as the group centre (BGG). In 
Section 3 we test the sensitivity of our caustic algorithm 
against the three group centre definitions. 



3.2 The mock group catalogue 

The mock catalogues were constructed by first populat- 
ing the dark matter halos of the Millennium Simulation 
( Springel et al.|2005 l with galaxies, the positions and prop- 
erties of which were predicted by the Bower et al. ( 2006 1 



description of the Durham semi-analytical model, GALFORM, 
and adjusted to match the GAMA survey luminosity func- 



tion of Loveday et al. (2011 1. To generate an all-sky survey 
to the GAMA depth, it is necessary to stack 11^ replicated 
copies of the simulation box to create a 'super-cube' (with 
sides of length 5500 h"^ Mpc). 

Following the prescription of Merson et al. (in prep) to 
place the halos and galaxies in the super-cube at the position 
at which they enter the light cone, galaxies are assigned 
their intrinsic properties, such as stellar mass and fiuxes in 
the appropriate bands. Galaxy properties are generally not 
interpolated between snapshots, with the exception of fiuxes. 
This is to avoid discontinuities at the snapshot boundaries. 
Having assigned the r-band apparent magnitudes, they are 
adjusted to perfectly match the redshift dependent GAMA 
galaxy luminosity function, following a simple abundance 
matching in the r-band. Finally the fiux limits of the GAMA 
survey are applied. We then take the all-sky mock catalogue 
and apply solid angle cuts to create the 9 mocks. 



4 CAUSTIC MASS ESTIMATES 

Before applying this algorithm to the actual group cata- 
logue, it is important to understand how well it performs 
when estimating masses for a set of mock catalogues that 
have been prepared alongside the G^Cvl. We therefore cal- 
ibrate our caustic mass estimation algorithm using a set of 
9 GAMA mock galaxy catalogues (described in section 2). 

In these mock catalogues, the true grouping of galaxies 
is known, so a well-informed calculation of their halo mass 
is possible. This acts as a benchmark for our caustic mass 
estimation algorithm and allow us to experiment with dif- 
ferent implementations of the method. |Serra et al.| ( |2011[ ) 
show that fine-tuning the parameters of the caustic mass 
algorithm {he, q, k) does not provide a considerable im- 
provement of its results. Instead we run the algorithm using 
different values for the group centre given in the G^Cvl: the 
centre of light, the iterative group centre and the bright- 
est cluster galaxy. This provides a useful way of testing the 
stability of the caustic algorithm to different definitions of 
the group centre, as well as confirming the conclusions from 
[Robotham et alT] ( |2011| ) of which definition is the most ap- 
propriate. 



4.1 Sensitivity to definition of group centre 



In Robotham et al. (20111, the iterative method seems to 



be particularly robust at picking out the group centre even 
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Figure 2. Example of a cut-out cone showing galaxies and groups with tab < 19.8 for the GAMA G12 field going from 174° ^ RA ^ 186° 
out to redshift z ~ 0.5. From left to right, the cones display all galaxies in this field, all groups found in G'^Cvl (coloured by their total 
group luminosity), all grouped galaxies and all ungrouped ones. 



in the presence of outliers, which tend to throw off the esti- 
mated group centre when using the centre of light method. 
The brightest cluster galaxy approach is also robust to out- 
liers, but analysis reveals that the iterative method recovers 
more group centres that match the mock groups. It precisely 
matches the group centre around 90% of the time, while the 
centre of light method is the poorest performer because it 
is the method that depends most explicitly on the group 
being recovered completely. For the caustic mass algorithm, 
we expect the most robust definition of the group centre to 
produce the most stable results; in other words, the vari- 
ance in the scale factors Ac for each multiplicity and red- 
shift bin should be minimal for the most stable group centre. 
This is confirmed when we run the caustic algorithm on the 
mock group catalogue and change only the location of the 
group centre from the iterative centre-of-light group cen- 
tre to the brightest group galaxy and then to the centre-of- 
light: cr^ (AitcrCon) = 0.192 whereas ^^(Abgg) = 0.201 and 
^"■^(^Col) ~ 0.361. We expect the tendency of the centre- 
of-light method to incorrectly define the group centre to be 
an outlying bright galaxy to throw off the placement of the 
caustics by deforming the density distribution of galaxies 
in redshift space. However, we also expect the caustic algo- 



rithm to be robust to minor perturbations in the placement 
of the galaxies in redshift space, as the difference in the vari- 
ances given above is minor. 

4.2 Sensitivity to definition of group extent 

Our second exploration involves changing the number of 
galaxies that we consider when calculating the mass for each 
group, i.e. artificially increasing group membership. We do 
this by extending the boundary of the group in redshift space 
to include some nearby galaxies. In redshift space, this ex- 
tension is defined as: 

\z — ^modl < Az; = 4 X max|z — Zmcdl (16) 

Spatially, we increase the maximum distance a given 
galaxy can be from the group centre. This effectively allows 
us to check whether the caustic algorithm is sensitive to 
other interloping groups. We find that including extra galax- 
ies that in some cases belong to other groups not only sys- 
tematically increases the mass estimates made by the caus- 
tic algorithm as one might expect, it also increases the mean 
spread of the results (defined as the ratio of the logarithm 
of the true and estimated mass) from (cr^) — 8.33 x 10^'^ for 
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Figure 3. Side-by-side comparison showing how the caustic mass estimation performs for a large mock galaxy group. In both panels, 
the black points represent the locations of galaxies in redshift space, the black lines show the contour fq{r, v) = k and the green 
lines are the caustics drawn along min{|»;u|, I"";!}- The left panel is a caustic fit only for galaxies present within this mock group and 
shows a clear example of the 'trumpet' distribution. The right panel includes nearby galaxies in redshift space out to it4iJioo and 
1^ — ^rnedl ^ 4 X max|2 — Zmedl whilst circling the original group members in red. In this case the trumpet distribution is lost and the 
caustics are artificially closed at the end of the sample. 



caustic mass estimates made with the groups as they are, 
to ((T^) = 0.0126. One expects the overall mass estimate to 
increase as a result of including more galaxies, but the in- 
crease in spread is unexpected: the caustic algorithm ought 
to work better with a greater number of galaxies. Instead, 
our algorithm is unable to correctly place the caustics in the 
redshift space diagram because our extended search cut in- 
cludes galaxies that are most likely inside other independent 
groups (Figure |5|. 

Visual examination of the redshift space plots in Figure 
[Sj of these extended group cuts demonstrates that the in- 
clusion of nearby galaxies disrupts the distinctive trumpet 
shaped distribution seen in ideal spherical groups. Figure [4] 
shows a mass comparison between the group catalogue data 
and the extended group cut. Shown are a subset of groups; 
those that with redshift between and 0.1, and with 5 to 
9 members (before including nearby galaxies). The caustic 
mass estimates made with the extended group cut (in dashed 
black lines) show a much greater spread, due to the reasons 
described above. 

This result highlights the importance of carefully deter- 
mining group membership when using the caustic algorithm, 
as the presence of galaxies that are not associated with the 
group being considered can have a catastrophic effect on the 
locations of the caustics, ultimately resulting in an incorrect 
mass estimate. By combining the caustic method with the 
friends-of-friends algorithm used to determine group mem- 
bership for the G'^Cvl and mocks we are able to signifi- 
cantly reduce the probability of the caustics being placed 
incorrectly due to the presence of interloper galaxies within 
a group. 



4.3 Application to mocks 

Using the original grouping from the G'^Cvl with the it- 
erative group centre, we calculate caustic mass estimates 
for every group in the mock catalogues and calculate Ac 
for a set of redshift and multiplicity bins in order to make 
these estimates median unbiased (Table 1). Figure [6] shows 
the results of this process for each multiplicity and redshift 
bin, comparing the distribution of caustic mass estimates (in 
black) to the dynamical mass estimates (red). Despite the 
caustic method performing best in very populated groups 
(A'^ ^ 200; of which there is only one such group in the 
G'^Cvl), demonstrated by the fact that the scatter of the 
distribution of mass estimates shown in Figure [6| decreases 
as a function of increasing multiplicity. We would normally 
expect the scaling factor Ac to always equal one. Instead we 
find that the caustic mass is systematically greater than the 
masses of the mock groups, i.e. Ac < 1. The only exception 
to this is for groups with two to four members, where the 
caustic algorithm is much more likely to fail to find appro- 
priate contours, and thus defaults to a specific value, arti- 
ficially adding a tail to the distribution of masses for that 
group subset. It is likely that due to the simulated galaxy 
groups in the mock catalogues not being perfectly spherical, 
the groups not being virialised (one of the basic assumptions 
of the spherical infall theory, on which the caustic algorithm 
is based, is that the group is virialised) the extra smoothing 
in our caustic introduced by using FFT to calculate the den- 
sity distribution fq{r,v), as well as the fact that the caustic 
method works best for galaxies with more than 200 mem- 
bers, Ac is not always equal to one as it should be in an ideal 
case. The overall variation in the scale factor Ac is roughly 
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a factor of 4.5, which is of the order of the range of scahng 
factors used to cahbrate the G'^Cvl dynamical mass esti- 
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Figure 4. The black dashed contours show the results of running 
the caustic algorithm on the mock catalogues and considering a 
larger number of galaxies around each cluster in the catalogues. 
By comparison, the solid black contours represent the mass esti- 
mates made when considering only galaxies known to be in each 
group from the simulation. Finally, the red contours represent the 
distribution of the dynamical mass estimates from the G^Cvl. All 
three distributions have been adjusted to the same median as the 
dynamical mass estimates, shown in the red contours. Each con- 
tour line contains 10, 50 and 90% of the points, and the dashed 
green lines are 2/5/10 times away from the median. 
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Figure 5. A projection of the mock group in Figure |3] showing 
all galaxies within the mock group circled in red, and all galaxies 
that are detected when the algorithm includes nearby galaxies out 
to it4/{ioo and \z — Zmedl ^ 4 X max|z — Zmedl The area shaded 
in blue represents the physical size of the simulation region at the 
median redshift of the group. At this range it is evident that the 
algorithm is including galaxies that are likely to belong to other 
separate groups. 



mates (see Table 3 in Robotham et al. (20111). In contrast 
to the dynamical mass estimate scaling factors, the caustic 
mass scaling factors vary less as a function of redshift, but 
are far susceptible to variations in group multiplicity. 

As expected, the least satisfactory results are for groups 
that are galaxy pairs, where there is little velocity informa- 
tion for the caustic algorithm to use. This explains the pres- 
ence of large tails in the 2 ^ A^poF ^ 4 panels of Figure [6] 
However, we must note that both methods tend to fail at 
these low multiplicities. For higher multiplicity and redshift 
cuts the caustic mass estimate PDFs shown in this figure 
agree extremely well with those of the dynamical mass esti- 
mates; indicating some correlation between both algorithms 
when it comes to performing badly for certain galaxy groups. 
This can be seen in the scatter plots as both contours tend 
to follow roughly the same profile, meaning that it is likely 
that a group that performs badly in one algorithm is likely 
to do so with the other. This is particularly visible for the 
10 ^ A^FoF ^ 19, 0.3 ^ zpoF ^ 0.5 bin where there is a 
secondary concentration of high mass groups that is present 
in both distributions. 

Figure [7] translates the information shown in Figure [6] 
into a set of density distributions where the ratio between 
the caustic and dynamical mass estimates and the known 
halo mass (shown in black and red respectively) is displayed 
alongside the ratio between both mass estimates (in blue). 
As seen on Figure[7j for groups with a mid-range multiplicity 
(5 ^ A^'foF 5; 19) the caustic mass estimates have a greater 
spread than the dynamical mass estimates; this is true across 
all redshift bins. For groups with A^foF 5; 4 the scatter is 
comparable across all redshifts. Both methods produce esti- 
mates that are within a factor of 2 in agreement with each 
other. The large tails seen in Figure [6] are visible here: the 
small 'bumps' forming on the right hand of the distribu- 
tion are seen in both distributions, once again demonstrat- 
ing that both methods tend to fail in similar ways. Despite 
the caustic algorithm being designed for high multiplicity 
groups, we are still able to make reasonable estimates of the 
group mass. 

Our aim is for the caustic mass estimate to recover the 
intrinsic halo mass of each group as accurately as possible. 
However, in practice and when applying the algorithm to 
the real data in the G^Cvl, we must run the algorithm not 
on intrinsic groups, but on groups defined by the friends-of- 
friends algorithm used in |Robotham et al. ( 2011[ ). An impor- 
tant test therefore is to see how the caustic mass estimation 
algorithm performs on bijectively matched groups drawn out 
from the mock catalogue using the same friends-of-friends 
algorithm, and the results of this can be seen in Figures [S] 
and [9] In both fi gures the distribution of data mimics that 
of the distributions for the caustic masses on the intrinsic 
mock groups, with these showing less scatter compared to 
the caustic masses of the FoF mock groups. Given that the 
design of the friends-of-friends algorithm used in |Robotham| 



et al. (2011 1 to construct the groups in both the mocks and 



the final group catalogue is to reject groups that have sig- 
nificant outliers, this greater agreement indicates that the 
caustic mass estimate may be slightly more sensitive to out- 
liers. 

One final check is to examine how the caustic mass es- 
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Table 1. Values for Ac for each group subset in both the tab < 19.4 and tab < 19.8 mock G'^Cvl mock catalogues using iterative 
group centres and with all galaxies in group. Including these numbers in Eg. |14| gives a median-unbiased estimate for the group mass. 



Mass estimates of true mock groups 
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Figure 6. Distribution of caustic (dynamical) masses for the intrinsic mock groups as a function of halo mass, drawn in black (red) for 
''AB ^ 19.4. For each panel, both mass estimates have been corrected to be median-unbiased. The contours represent areas containing 
10, 50 and 90% of the groups and the green lines are regions where the mass estimate is 2/5/10 times off the true mass. Of particular 
interest is the tendency for both distributions to follow each other very closely, particularly when over or under estimating the true halo 
mass. 
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Figure 7. Probability distribution functions of median-unbiased values of log 



for the rA.B < 19-4 sample, with caustic masses 



drawn in black, and dynamical masses in red. The blue line is the PDF of log . , " and highlights the agreement between the two 
methods. The green dashed lines indicate regions that are factors of 2/5/10 away from the 'true' mass. The difference between scatter in 
the caustic and dynamical mass estimates is often minimal, with the caustic method showing less scatter for groups with 5 ^ A''foF ^5 19. 



timates perform as a function of the quality of the grouping, 



as defined in Robotham et aL (2011 1. The total quality pa- 



rameter, QTot, is a measure of how significantly matched 
individual groups are. The best two-way matching group is 
defined as being the one which has the largest product for 
the relative membership between the recovered FoF group 
and the closest matching mock group (see Eqs. (12) to (14) 
in Robotham et al. (2011 \) . In Figure 10 we show how the 
two mass estimates behave as a function of Qto*. From the 



left panel we are able to see that the caustic mass estimates 
behave in just the same way as the FoF mock dynamical 
mass estimates, performing well after a total quality factor 
of about 0.2. This is further demonstrated in the right panel 
of this figure, where the ratio between the caustic and dy- 
namical mass estimates for the FoF mock groups remains 
close to unity as a function of Qto*, with a small tendency 
for the caustic mass to be systematically greater than the 
dynamical mass as Qtoi approaches 1. 
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Figure 8. As Fig. [6] but for bijectively matched groups identified by the friends-of- friends algorithm instead of the true known intrinsic 
grouping. The coloured lines represent regions containing 10, 50 and 90 % of groups. The black contours compare the caustic mass to 
the FoF halo mass, and FoF dynamical mass estimate to the FoF halo mass estimate are shown by the red contours. 



Based on these explorations, we chose to create our 
caustic mass estimates for the mock groups in the catalogue 
using group centres obtained with the iterative CoL rejection 
method, and to only include galaxies considered to belong 
to each group. To generate each set of results, we select the 
same redshift and multiplicity cuts used in |Robotham et al.| 
( |2011[ ) We calculate the scaling factors for each bin as 
the necessary value to ensure that the median of the ratio 
of the mass estimate to the true mass is unity. Values for Ac 
are listed in Table [T] 



4.4 Application to G^Cvl 

With the appropriate scale factors in Table [l] we run the 
algorithm on the actual group catalogue itself, using the 
same redshift and multiplicity bins, the same group centre 
and only galaxies within each group. This provides us with 
a full set of halo mass measurements for every group in the 
catalogue that are complementary of those already existing 
in the G^Cvl. 

Both the caustic and the dynamical mass estimates are 
adjusted by their appropriate scale factors, and then com- 
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Figure 9. As Fig. [7] but for the bijectively matciied FoF groups from tlie mock galaxy catalogues. The different lines show median- 
unbiased values of log for the tab < 19-4 sample. The black line shows the ratio between the caustic mass and the FoF halo mass, 
while the red line shows the ratio between the FoF dynamical and halo masses. As before, the blue line shows the distribution of the 
ratio between the caustic mass and the FoF dynamical mass. 



pared. The results of this can be seen in Figure [TT] for the 
i^AB < 19.4 mag limited sample. By calculating the ratio of 
the caustic to the dynamical mass, we show that 90.8% ± 
6.1% groups have a caustic mass estimate that is within a 
factor of two of the dynamical mass estimate. In a vast ma- 
jority of cases the two mass estimates agree very well with 
each other, particularly for groups at high redshift (bot- 
tom row). In all cases the scatter of the data is minimal, 
with the scatter decreasing with increasing group member- 
ship through a combination of better quality mass estimates 



and the smaller number of galaxies present in these sub- 
sets. When interpreting these results, one must bear in mind 
that the caustic algorithm has been designed with popu- 
lous groups in mind; ideally with more than 200 members. 
The G'^Cvl has exactly one group that fits this criteria, at 
264 members. Of 12,200 groups in the G^Cvl, 10,813 of 
these have between 2 and 4 members and there are only 68 
groups that contain more than 20 galaxies. We also observe 
a disparity between the number of compact groups in the 
mocks compared to those in GAMA; this is further discussed 
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Figure 10. Left: A comparison of the performance of the two mass estimation methods as a function of the grouping quaUty parameter, 
Qxot- Black points correspond to caustic masses, and red ones to dynamical masses. Both methods perform worst when Qxot is close to 
0, but quickly recover at approximately 0.2. The green lines show a rolling median in each 0.1 bin for the caustic mass (solid line) and 
for the dynamical mass (dashed line). Right: The ratio between the caustic and the dynamical mass estimates as a function of Qxoti 
with the rolling median in each 0.1 bin overplotted in green. 



in Robotham et al. (20111. One also expects small-number 



statistics to affect the results for these larger groups: this is 
particularly evident in Table [2] which shows the fraction of 
caustic mass estimates within certain factors of the dynami- 
cal mass estimate. As group multiplicity increases these frac- 
tions become less reliable, as the number of groups present 
in each bin drops sharply. 

In order to make a comparison between the groups and 
mocks as fair as possible, we must account for the fact that 
velocity uncertainties have not been included in the mocks 



as this directly affects the kernel, through Eqn. (111. This 
issue is addressed in [Robotham et al.| ( |2011| ), whereby all 
groups with ^ 130 km s~ are removed from any com- 
parisons between the groups and mocks, arguing that below 
this cut the velocity dispersion of a group would be signifi- 
cantly affected by this uncertainty and that the dynamical 
mass estimate is directly proportional to a^. Based on this 
result, we apply the same velocity cut to the group and mock 
catalogues. Figures [Tl] and [12] display the resulting compari- 
son between the ratio Mc/Mdyn of the caustic and dynamical 
masses for the groups (black) and mocks (black). There is 
still some difference between the scatter of the mass esti- 
mates for the groups and the mocks and the disagreement 
is most noticeable in cases where the caustic mass is greater 
than the dynamical mass. However, if we do not include 
this velocity cut and include the mock groups with poorly 
defined masses, the discrepancy between the two distribu- 
tions increases, with Mc/Mdyn for the groups showing much 
less scatter than for the mocks. We note that when calcu- 
lating the mass estimates for the bijectivel y m atched mock 
groups with 130 km'^ s~'^ for Figure 

and ar to tend to 0, as Eqn. [9] diverges if hi 
hopt = 0. In much the same way as done in Table |2] we are 
able to calculate that the same percentage of caustic and 



we allow av 
when 



dynamical mass estimates fall within 50% of the FoF halo 
mass, implying that both algorithms have correctly been 
calibrated against the mocks. When comparing the ratios of 
both mass estimates to the known halo mass for bijectively 
matched groups, we find that these ratios are virtually the 
same. From this we can conclude that any final discrepancies 
in the mass estimates shown in Figures [TT] and [T2| must be 
down to an intrinsic difference between the mocks and the 
real data, and that both methods are properly calibrated to 
output results that are median unbiased with respect to the 
halo masses. 

Taking these considerations into account, these new 
mass estimates appear to agree with the existing dynam- 
ical mass estimates for the G^Cvl. While both methods 
utilise the same velocity and positional information in dif- 
ferent forms, the dynamical mass method is not as sensitive 
to the full 2D velocity profile of each group. This is caused 
by the dependence of v and r to each other, which is another 
effect that is considered in the caustic method, whereas the 
velocity dispersion and group radius used to calculate Mdyn 
are treated independently to each other. 



5 CONCLUSION 

Using the caustic mass estimation algorithm introduced by 
DG97 and D99 we have provided complementary caustic 
mass estimates for the G^Cvl. We have calibrated our im- 
plementation of the algorithm by running it on mock GAMA 
group catalogues where the intrinsic grouping is known; and 
for bijectively matched friends-of-friends mock groups. This 
allows us to derive a scaling factor Ac for bins of redshift and 
multiplicity which provides us with median-unbiased mock 
mass estimates, and we find that in all cases there is very 
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Figure 11. Distribution of the caustic masses compared to dynamical masses for real galaxy groups from G'^Cvl. The dashed green lines 
are as Figure [6] and the solid black lines represent regions containing 10, 50 and 90% of the data. There is a good correlation between 
the caustic and dynamical mass estimates as already shown in Fig [8] and [9] for the mocks. 



good agreement between the caustic mass and dynamical 
mass estimates. Of great note is the tendency for both dis- 
tributions shown in Figure |6] to match each other extremely 
well, implying that both algorithms perform equally well and 
equally badly for mock groups: that both estimate methods 
fail for low multiplicity {2 ^ N ^ 4) groups highlights an 
interesting systematic bias present in both methods that 
warrants further study. 

Having calibrated the algorithm on the mock cata- 
logues, we apply it to the G3Cvl and obtain caustic mass 
estimates for real groups as shown in Figure [TT| As with the 



mocks, we demonstrate that both mass estimates are gen- 
erally consistent with each other: on average, 90.8 ± 6.1 % 
of groups have caustic and dynamical mass estimates that 
agree to within a factor of two. This is strong evidence for 
the reliability of the caustic mass estimation method, par- 
ticularly when considering the less than ideal conditions the 
algorithm faces when working with so many groups of low 
multiplicity. Despite being designed to work with groups 
with over 200 members, the caustic method is able to suc- 
cessfully determine accurate mass estimates for groups with 
as few as two members, highlighting the adaptability and 
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Figure 12. As in Figure [7] this figure shows the PDF of the logarithm of the mass ratio between the caustic and dynamical mass 
estimates for all groups with cr^ ^ 130 km^ s~^. The blue line shows the same distribution but for the bijectively matched FoF mock 
catalogues with the same velocity cut to illustrate the smaller scatter between the two methods for the real data. 



strength of this method when combined with a friends-of- 
friends algorithm. We stress the importance of applying the 
caustic mass algorithm alongside a friends-of-friends algo- 
rithm; by accurately determining the membership of any 
given group we reduce the chances that the caustics may 
be located incorrectly. The tightness of the contours shown 
in the first column of Figure [ll] and the large fraction of 
accurately determined masses given in Table [2] are very con- 
vincing demonstrations that the caustic method is able to 
perform well even for very low multiplicity groups. 

We wish to draw particular attention to the result 



shown in Figure [12] that the ratio between the caustic and 
dynamical mass estimates in the groups is very sensitive to 
the application of a velocity cut at ^ 130 km^ s~^, as re- 
quired by the fact that the mocks (and gene the mock group 
catalogues) have not been analysed with realistic velocity er- 
rors. By discarding groups whose velocity dispersions would 
be badly affected by the inclusion of velocity errors, we en- 
sure that we apply a fair comparison to the group catalogue, 
which does include velocity errors (and these velocity errors 
play a crucial part in the calculation of the kernel on which 
the caustic algorithm depends). Once this cut is applied. 
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2 s£ AfpoF ^4 5 ^ ATpoF ^9 10 s£ TVpoP ^ 19 20 s£ TVpoP ^ 1000 





2 


5 


eiN2) 


2 


5 


e{N2) 


2 


5 


4^2) 


2 


5 




^ 2PoF ^ 0.1 


0.878 


0.970 


0.010 


0.862 


0.983 


0.030 


0.958 


1.0 


0.030 


0.769 


1.0 


0.152 


0.1 < ^FoP 0.2 


0.898 


0.977 


0.005 


0.836 


0.904 


0.020 


0.949 


1.0 


0.020 


0.978 


1.0 


0.022 


0.2 ^ 2FoP s£ 0.3 


0.904 


0.979 


0.006 


0.794 


0.997 


0.028 


0.952 


1.0 


0.028 


1.0 


1.0 





0.3 s£ zpoP s£ 0.5 


0.909 


0.978 


0.009 


0.904 


0.979 


0.034 


0.889 


1.0 


0.118 


1.0 


1.0 






Table 2. Fraction of groups that have a caustic mass estimate that is within a factor of 2 and 5 from the dynamical mass estimate. The 
error e(N2) is defined as yJl/N2 — 1/Ntot where N2 is the number of groups within a factor 2 and A^tot is the total number of groups 
for a given bin. The number of groups present in each multiplicity bin drops sharply (from thousands of groups to 10 or less) after the 
first multiplicity bin, drastically lowering the relevance of these statistics. 



we show that the distribution of ratio between the caustic 
and dynamical masses is consistent between the mocks and 
groups, particularly for groups with 2 to 4 members. 

These caustic mass estimates serve to add further cred- 
ibility to the dynamical mass estimates, as the caustic al- 
gorithm is sensitive to the escape velocity profile of a given 
group out to its full extent. That both dynamical mass and 
caustic mass estimates for the mocks and the groups give 
such similar results despite their relative independence from 
each other is a positive result that not only reinforces the 
dynamical mass estimates, but also demonstrates the power 
of the caustic approach when combined with a friends-of- 
friends algorithm. 

The G^Cvl will be made publicly available on the 
GAMA website as and when the associated redshift data 
is made available. Those wishing to use the group cata- 
logue before this time should contact ASGR at asgrQst- 
and.ac.uk. The caustic mass estimation algorithm devel- 



oped by the author has been written in R (R Development 
Core Team|20TT| . A non-GAMA specific version of the algo- 
rithm is freely available to use and can be downloaded from 



http: / /www. gama-survey.org/ pubs / 
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